On the Theory of Quantum Oscillations of the Elastic Moduli in Layered Conductors 
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In this paper we study theoretically how the local geometry of the Fermi surface (FS) of a layered 
conductor can affect quantum oscillations in the thermodynamic observables. We introduce a con- 
crete model of the FS of a layered conductor. The model permits us to analyze the characteristic 
features of quantum oscillatory phenomena in these materials which occur due to local anomalies of 
the Gaussian curvature of the FS. Our analysis takes into account strong interaction among quasi- 
particles and we study the effect of this interaction within the framework of Fermi-liquid theory. 
We show that singularities in the density of states of quasiparticles (DOS) on the FS which occur 
at low temperatures near the peaks of its oscillations in a strong magnetic field can be significantly 
strengthened when the FS of the layered conductor is locally flattened. This can lead to magnetic 
and lattice instabilities of a special kind which are considered in the final part of the work. 

PACS numbers: 71.10 Pm, 73.40 Hm, 73.20 Dx 
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I. INTRODUCTION 

Many superconducting materials with large critical pa- 
rameters created in the last two decades are layered struc- 
tures with metallic-type conductivity (e.g organic metals 
of the cv-(BEDT-TTF) 2 MHg(SCH) 4 group). A charac- 
teristic feature of these materials is a strong anisotropy 
of the conductivity in the nonsuperconducting state: the 
conductivity in the layer plane is much higher than in the 
direction normal to the layers. It is common to assume 
that an anisotropy in electrical conductivity is a manifes- 
tation of the quasi-two-dimensional nature of the energy 
spectrum of the charge carriers in layered conductors. It 
follows from experiments that the Fermi surface (FS) of 
such conductors is a system of weakly rippled cylinders 
(isolated or connected by links) whose axes are perpen- 
dicular to the layers. The experimental data (see, e.g. 

Refs. [1EHSIEIS00I3) su PP° rt this assumption. 

The FS of a conductor with a quasi-two-dimensional 
energy spectrum can be described by the following equa- 
tion 



k 

E(p) = Ek (P* ' Py ) C0S 



(1) 



where p is the electron quasimomentum; Ek(p x ,Py) are 
coefficients with dimensions of energy, p z is the projection 
of the quasimomentum on the direction normal to the 
layers, and a is the distance between the layers. If we 
ignore the anisotropy of the energy spectrum in the layer 
plane, instead of (1) we can write the simpler equation: 



E p = 



— + 5>(^)cos — 

fc=l 



(2) 



where p_L is the projection of the quasimomentum on the 
layer plane, and m± is the effective mass corresponding to 
the motion of the quasiparticles in that plane. Equation 
(2) describes an axially symmetric open FS with the axis 



of symmetry (z axis of the chosen coordinate system) 
directed along a normal to the layers. The profile of the 
longitudinal section of the FS corresponding to Eq. (2) 
is a periodic function of period 2po = ^irh/a. 

Experimental data concerning quantum oscillations in 
various characteristics of metals were widely used as 
the instruments of study of their electron energy spec- 
tra. Quantum oscillatory phenomena were repeatedly ob- 
served in organic layered conductors (Refs. yj, y, an d 
|lfj|). Quantum oscillations in thermodynamic observ- 
ables in these materials can exhibit some characteristic 
properties which occur due to the quasi-2D character of 
the energy spectrum of charge carriers. For instance it 
was shown in Ref. that the magnetic susceptibility of 
the noninteracting 2D electron system with open FS can 
exhibit sharp maxima due to the orbital effect. However, 
a systematic theory of quantum oscillatory phenomena 
in organic metals which takes into account the quasi-2D 
character of their electron spectra and a strong correla- 
tion among the electrons is not well developed at present. 
In this paper we give such theoretical analysis and we ap- 
ply our results to study the quantum oscillations of the 
thermodynamic characteristics, especially of the elastic 
moduli. 

At first, we derive the expressions for the electron con- 
tributions to the elastic constants corresponding to the 
compression in the direction perpendicular to layers (Ao) 
and to the compression applied along the layers (A) in 
the presence of external magnetic field B. In the geome- 
try where the "z" axis is directed along the normal to the 
layers, and assuming that the system is axially symmetric 
our elastic constants Ao and A correspond to the electron 
contribution to the elastic moduli C33 and en = C22 in 
Voigt notation. A possible instability of the lattice arises 
near the peaks of the oscillations of DOS in a quantizing 
magnetic field. Therefore we analyze these oscillations of 
DOS in detail, and we show that they can be strongly 
affected due to the interaction among the carriers. In 
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the following section of the work we introduce a concrete 
model of the FS which permits us to analyze how the 
local geometry of the FS can strengthen the singulari- 
ties in DOS which can give rise to the considered lattice 
instabilities. 



2. MAIN EQUATIONS 

To simplify further calculations we do not consider the 
deformation interaction between the electrons and the 
lattice. Using the generally accepted assumption the ef- 
fect of the electrons on the lattice arises due to a self- 
consistent electric field, which occurs under deformation. 
Also, the deformation of the lattice causes the occurence 
of an additional non- uniform magnetic field b(r). The 
presence of these fields leads to a redistribution of the 
electron density N. The local change in the electronic 
density 6N(r) equals: 



6N(r) 



■^(r) + — b(r) 



-N? 



dB 



b(r) 



(3) 



Here, e is the absolute value of the electron charge, B is 
the external magnetic field. 

The magnetic field b(r) satisfies the equation: 

curl b(r) = 47r curl £M(r) 



47r curl 
0. 



dM 



eip(r) 



dM 
dB 



b(r) 



div b(r) 



(4) 



Here, M is the magnetization vector; £ is the chemical 
potential of charge carriers; (p(r) is the potential of the 
electric field, arising due to the deformation; N$ is the 
"bare" density of states (DOS) of quasiparticles on the 
Fermi surface of the layered conductor: 

dU 
dE v 



V — 



(5) 



f v is the Fermi distribution function for the quasiparticles 
with energies E v . 

Within the framework of quantum phcnomenological 
Fermi-liquid theory the quantity N£ is the DOS of the 
quasiparticles renormalized due to the Fermi- liquid inter- 
action among them (see Ref. 0]): 



n; 



i E v — E,,i 



< W '(-q)^'f(q)|g=o- (6) 



and the renormalized operator of electron density 
n* vv , (— q) is connected with the "bare" operator n vv i{— q) 
by the relation: 



Ku'i-l) = »w(-q) + 



E„ 



E u 



•FT. 



-q) 
(7) 



where F^f 2 are the matrix elements of the Fermi-liquid 
kernel [See Eq. (13) below]. 

The relations (3), (4) have to be complemented by the 
condition of electrical neutrality of the system: 



SN(r) +eNdiv u(r) = 0, 



(8) 



where u(r) is the lattice displacement vector, TV is the 
electron density. It follows from equations (3), (4) and 
(8) that: 

f dM N 1 

curl{(l - 4vrx)b(r)} = -47rcurW -^-^rdiv u(r) \ . 

(9) 

The set of simultaneous equations (3), (4) and (8) was 
first presented in Refs. |l4j and [15l |. It allows us to 
exclude b(r) and to express the potential <p(r) in terms 
of the lattice displacement vector. As a result we can 
derive the expression for the electron force F(r) acting 
upon the lattice under its displacement by the vector 
u(r). For the conductor with axially symmetric FS in 
the magnetic field B directed along the symmetry axis, 
the force F(r) equals: 

F(r) = A b (b V(Vu(r))) + A[b X [V(Vu(r)) x b ]]. 

(10) 

Here bo is the unit vector directed along B. We assume 
that the magnetic field B is directed along the normal to 
the layers. Then bo coincides with the normal vector to 
the layers n. It follows from Eq. (10), that under longitu- 
dinal (parallel to bo or perpendicular to layers) direction 
of gradients of the deformation tensor the electron con- 
tribution to the corresponding elastic modulus is equal 
to Ao, and under their direction across bo (in the plane 
of layers) it equals A. The elastic constants Ao and A are 
described by the expressions: 



An 



N 



A = Ao 1 



1 - 4tTX|| 



(11) 

(12) 



Here, x\\ is t ne longitudinal part of the magnetic suscep- 
tibility; xc — {9M z /d()(dC/9B). More careful treatment 
carried out before including the deformation interaction 
(see e.g. Ref. [16|), do not give a qualitative change of 
results. If we will take into account such interaction, it 
will lead us to a replacement of the electron density N 
in the expression (11) for Ao by the quantity N(l + a), 
where the value of the dimensionless constant a depends 
on the corresponding component of the tensor of defor- 
mation potential averaged over the FS. The constant N 
in Eq. (12) also has to be replaced by N(l + a*) where 
a* is a dimensionless constant depending on the averaged 
over the FS deformation potential. 

As follows from Eq. (11), in our case, when we neglect 
the deformation interaction the quantity Aq coincides 
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with the compression modulus of the electron liquid. The 
structure of the quantity A is more complicated. Besides 
the contribution, connected with the electron compres- 
sion itself, A also contains the contribution arising due 
to the magnetostriction (the second term in the expres- 
sion enclosed in the brackets). It is shown below that at 
low temperatures when 6 <C 1 {0 = 2tt 2 T '/ 'HQ; T is the 
temperature expressed in energy units, Q is the cyclotron 
frequency of the charge carriers) there is a significant dis- 
tinction between the oscillating corrections to the elastic 
moduli Ao and A arising due to the magnetostriction. In 
the considered case (B||n) the velocity of the longitudi- 
nal sound, propagating along the normal to the layers is 
proportional to the square root of Ao, and when it prop- 
agates in the plane of layers - to the square root of A. 
Therefore, the distinction between Ao and A due to the 
magnetostriction will be displayed in the quantum oscil- 
lations of the velocity of sound. Above the low temper- 
ature range {9 ~ 1) the contribution connected with the 
magnetostriction becomes negligible and the distinction 
between Aq and A is smoothed. 



LOW TEMPERATURE QUANTUM 
OSCILLATIONS OF THE DOS 



Here, we calculate the oscillating corrections to the 
renormalized DOS Nt which is defined by Eq. (6). The 
calculation is carried out for an axially symmetric FS of 
an arbitrary shape and this needs a special consideration 
because the phenomenological Fermi-liquid theory is well 
developed only for isotropic model of a metal. We start 
from the well known expressions for matrix elements of 
the Fermi-liquid kernel F^Jf 2 |l7j : 



pVlV2 _ 

vv' 



w aia2 S ,S 4- 4?/;" lQ2 fs /s 



1717' a (7\l72 ) 1 



(13) 



where a is the set of the orbital quantum numbers, a is 
the spin number, s is the spin operator. 

The FS under consideration is axialy symmetric. 
Therefore Fermi-liquid functions ip(p, p') and ip(p, p') on 
the FS depend only on the cosine of the angle 9 between 
the vectors pj_ and p' ± and on the projections of the 
quasimomenta on the symmetry axis p z ,p z . We can rep- 
resent the function </?(p, p') as a sum of even and odd 
contributions in cos 9 : 

<p(p,p') = fo{Pz,p' z ,cos9) + cos9tpi( Pz ,p' z ,cos9), (14) 

where ipo (p z , p' z , cos 9) and tp\ (p z , p' z , cos 9) are even func- 
tion in cos#. By virtue of the FS invariance under re- 
placements p — > p and p' — > — p', the functions 
ipo(p z ,p' z ,cos9) and ipi(p z ,p' z ,cos9) should not vary un- 
der simultaneous change of sign of p z and p' z . It enables 
us to separate each of the mentioned functions into the 
parts which are even and odd in p z p' z and to rewrite Eq. 
(14) in the form: 

v(P,P') = f oo + PzP' z f oi + cos6>(</Jio +p z p' z (pn). (15) 



Here, the functions ipoo, tpm> ¥>io and ip\\ depend onp z ,p' z 
and cos# and are even functions in all their arguments. 
We can write a similar expression for ijj(p,p'). 

To proceed we expand functions included into Eq. (15) 
in Fourier series in the angles <E> and The functions 
are even in cos 9, therefore the expansions have the form: 

oo 

Vlk ( Pz ,p' z ,cos9)= ^\pz,p' z )e 2ls{ ^' ] . (16) 

s— — oo 

Here i,k = 0,1. Coefficients ^j^' (p z ,p z ), in the expan- 
sions (16) are even functions in both arguments. Be- 
sides, they should not change when the arguments are 
interchanged. We can simplify further consideration as- 
suming that: 



1 



9o 



(17) 



where both factors are even functions, and g$ is the 
"bare" DOS of the quasiparticles on the FS in the ab- 
sence of the magnetic field B. 

The expansions of the functions Tp£. (p z ,p' z , cos 6) are 
similar to Eq. (16). The coefficients in these expansions 
also can be presented as: 

i>t s \p*,Pz) = ^Ql s{ (Pz)Q? k sl (p'z), (is) 

9o 

To calculate the quantity N£ under the conditions of 
the semiclassical quantization (7 3> 1 where 7 2 = 
we can replace the matrix elements p^? 2 and ^aa" 2 by 
their semiclassical analogs. The analog for f^f 2 is the 
Fourier component in the expansion of the function: 



p(q;p,p') = ex P 

x exp 



icq y pA 
eB 

icq v p' A 
eB 



sin $ 



sin $' 



¥>(P,P') 



(19) 



in the Fourier series in the angles $ and <!>'. Here, <p(p, p') 
is the semiclassical Fermi-liquid function. With the help 
of the known identity for the Bessel functions: 

00 

exp(±iz sin 

m— — 00 

we arrive at the result: 
1 



W(q;p,p') 



\ " 



T ( cq y p± 

m— — oo 



J m > ( £ |f= L )j[ ^cxp[z(m-;]$] 

/•Z7T 

/ d$'exp[i(m'-Z')$V(p,p')- (21) 
Jo 



rn' — — 00 
/•2tt 
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Substituting the expression (14) for (p(p,p') into Eq. 
(21) and taking into account the expansions (16) we ar- 
rive at the final result: 



¥>H'(q;P,p') 

1 oo 



S = — OO 



v ' j ( cq y P± \ T ( cq y p' ± \ 

, x f c 1yP± \ T ( c( lyP'i 

+ ^-2s+l I J J Z'-2s+l 



V eB 



(22) 



In the expansions of ^/(q; p, p') the functions -^'^'(.Pz) 
and -Pj? 8 ' (p' 2 ) in (22) have to be replaced by the functions 



Qll s \ Pz ) and Q[ z k s \p' z ). The derived expressions (22) en- 
able us to study any Fermi-liquid effects for conductors 
whose FSs are axially symmetric. 

Here, we consider the most interesting case when the 
gradients of the deformation of the lattice are directed 
across bo(q_Lb ). Introducing the notations 2N — 



2s 



2N' = ni 



n 2 ; 



I = n — n': V = n\ — n<i 



we can write the expansions for the matrix elements 
(p(q;Nlp z ; N'l'p' z ) and V(q; Nlp z ; N'Vp' z ). We arrive at 
these expansions replacing the quantities p±,p'^ in the 
arguments of Bessel functions in Eq. (22) by their quan- 
tum analogs hy/2N + 1 and hy / 2N' + 1. 

The functions Pj£ s ' (p z ) and (p z ) are even, and the 
matrix elements n vv r(—q) are diagonal in spin number a. 
Owing to this, one can take into account only the first 
term in the first square brackets in Eq. (22) in calcu- 
lation of the quantity N£. The contributions from the 
other terms will be equal to zero. Thus in the following 
calculations we can assume that ip(q;Nlp z ; N'l'p' z ) and 
Nlp z ; N'l'p' z ) have the form 



1 



<pfam PM] N'i' P ' z ) = - 



P^{P*)P^{P' Z 



90 .±rt. 



|2s|, 



<P(q]Nlp z ; N'l'p' z 



x J ( _ 

x J v 
1 

.9o 

X Ji_ 2 
x Jl>- 



2 s 



-2 s 



eB 

^V2WT 

eB 



E 

S — — oo 



Q&''(p»)Q|£V,) 

hcq 



eB 
hcq 
~e~B 



V2N + 1 



V2N> 



(23) 



Further we omit the lower indices of the functions 



p al^(Pz) and QocTfoO for simplicity. 

Under conditions of the semiclassical quantization, 
when 7 > 1, we have for the arbitrary single particle 
operator $,,„/ (see Ref. 0): 



,128b 



E 



E,, — E,,i 



■EE/* 



m ± (p z W (p z ) 



+ 5 

Iq Qaa' 



(Pm) 



.(24) 



Here, m±(p z ) is the cyclotron mass which is assumed 
to be a constant in further calculations; p m is the value 
of the longitudinal component of the quasimomcntum, 
corresponding to the mth external cross section of the 
FS. The functions A and A' describe oscillating correc- 
tions which appear in quantizing magnetic fields. Their 
form depends on the particular character of the energy- 
momentum relation for the quasiparticles. We will ana- 
lyze these quantum oscillating corrections below within 
the framework of our model for the electron energy spec- 
trum in organic metals. Integration with respect to p z in 
Eq. (24) is carried out within the limits determined by 
the shape and size of the FS. 

Using the asymptotic formula (24), and the expansions 
(23) we can obtain the expression for the renormalized 
DOS: 

(1 - a ) 2 [A + /3(A 2 - A' 2 )] 



N? = go 



l-a 



1 + {a + /3)A + a/3(A 2 - A' 2 ) 



(25) 



Here: 



o-o 



Ao = 
A = 
A = 
B = 

An = 



A 



l + Ao' 
A* n 



1 o 



1 + A 



do 



(3 = 



A Q 



B* 



o 



1 + Bn 



27T 2 h 6 go 
P°(0)mj 
2^h 3 g a 

2^¥g~ 

m± 

2TT 2 h 3 g 
[P»] 2 ; 



P°(pz)d P: 

P°(p z )d Pz , 

P\p z )fd Pz . 

Q°(p z )} 2 d Pz ; 
B* = [Q»] 5 



(26) 

Neglecting terms proportional to the differencies A 2 — A' 2 
in the numerator and the denominator of the Eq. (26) we 
can simplify the expression for the renormalized density 
of states: 



9o 



1 - a + 



(l-5o) 2 A 
l + (a + /3)A 



(27) 
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The expression (27) is easily generalized to cover the case 
when the Fermi surface has several extremal cross sec- 
tions. In this case we obtain 

£(l-5 m ) 2 A m " 

= 5o 1 - «o + 1 i a u = 9a(l-a +K). 

1 + LK + Pm)A m 

m 

(28) 

Here, summation over m is carried out over the extremal 
cross sections: 



A* 



A, 



B* 



(3 m 

l + A Q ' "" 1 + A ' 1 + B a 

A* m = [P°( Pm )} 2 ; B* m = {Q°( Pm )} 2 ; 

P u { Pm ) I P°( Pz )d Pz . 



A 



m ± DO, 



2ir 2 h 3 



(29) 



Neglecting all terms arising due to the interaction among 
quasiparticles, we arrive at the result for the "bare" DOS: 



N ( = g (l + A). 



(30) 



The oscillating term A describes quantum oscillations 
of the "bare" DOS. This term is small compared to unity 
when the temperature is moderately low (9 ~ 1). When 
9 « 1 the magnitude of the function A can reach val- 
ues of the order of unity. Comparison of Eqs. (28) and 
(30) shows that there are pronounced distinctions be- 
tween the "bare" DOS N c and the DOS N£ renormal- 
ized due to the interaction among quasiparticles within 
low temperature range. Under 9 <C 1 the amplitude and 
the form of the function K which describes the quantum 
oscillations of N£ can differ significantly from the corre- 
sponding characteristics of the function A. In particular, 
the denominator in Eq. (22) can turn zero at the peaks 
of the oscillations. Thus, the function K has singularities 
caused by the Fermi-liquid interaction. 

The specific features of the low temperature quantum 
oscillations in renormalized density of states have to be 
manifested in oscillations of the thermodynamic observ- 
ables. We consider, for example, the oscillating contribu- 
tion to the magnetic susceptibility for a uniform field B. 
To simplify the following calculations we assume that the 
unique extreme cross-section of the axisymmetric FS is 
located at Pz — 0. We remark here that possible singular- 
ities in the longitudinal magnetic susceptibility near the 
peaks of quantum oscillations of the renormalized DOS 
were first analyzed in Ref. f° r an isotropic electron 
liquid (see also Re f. Ilfll l. Generalizing the correspond- 
ing result of Ref. [18j, we can show that in the case of 
the axisymmetric FS, \\\ is proportional to the oscillating 
function K' : 



K' 



K 



l + [l-a (<S)]K 

A(0) 




l + [l + ao(0)-ab(0)+A>(0)]A(0)' 



(31) 



B x B B 2 

FIG. 1: Dependence of the magnetic susceptibility \\\ on 
the magnetic field at T — 0. The domain of the values of the 
magnetic field where 4n\ > 1 is section-lined. 



Within the low temperature range the amplitude of the 
oscillations in A(0) might increase to such an extend that 
it reaches the values of the order of unity. Under these 
conditions the form of the de Haas-van Alphen oscilla- 
tions is determined by the denominator in Eq. (31) to a 
great extent, and it depends critically on the Fermi- liquid 
interaction. When [1 + a (0) - a (0) + /3 (0)] < the 
denominator in Eq. (31) can go to zero at the peaks of 
the oscillations and the susceptibility x\\ would increase 
beyond all bounds. It would lead to the magnetic insta- 
bility (x|| < l/47r) at the peaks of the oscillations (see 
Fig. 1). Thus, the low temperature analog for the Shoen- 
berg effect can be observed. It follows from the results 
of Ref. [18j as well as from our expression (31) that the 
magnetic susceptibility \\\ of noninteracting 3D electron 
gas can exhibit a series of maxima of a finite magnitude 
when temperature goes to zero. G. Montambaux et al. 

came to a similar conclusion irrespective of the ear- 
lier results of Ref. . 

These singularities of the renormalized DOS N£ orig- 
inating due to the Fermi-liquid interaction can give rise 
to specific features in the quantum oscillations of ther- 
modynamic observables in layered conductors. However, 
their manifestations would depend on the local geome- 
try of the FS near its extreme cross sections. To analyze 
these anomalies in more detail we introduce a concrete 
model of the FS for layered organic metals which allows 
us to consider a broad class of the FSs of various profiles. 



4. MODEL AND RESULTS 

The usual approach in theoretical studies of electron 
properties of layered conductors is to keep only the first 
few terms in the sum over k in Eq. (2). As a rule (see 
Refs. [jfl EH E2T | , only the first term is taken into ac- 
count. This corresponds to results obtained in the tight- 
binding approximation. Here, we use a different approach 
to describe the electron energy spectrum of the charge 
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carriers in layered conductors, whose FS is denned by 
the equation 



E(p) = -±- - T]v p E — 
2m ± \po 



(32) 



where vo = (2(/mi) 1 / 2 , E(p z /po) is an even function 
periodic in its argument p z /po with a period equal to 
2, and i] is a dimensionless parameter characterizing the 
rate of rippling of the FS. The quantity — r]VopoE(p z /Po) 
is the sum of the trigonometric series in (2) . By selecting 
the type of this function we can obtain FS's shaped as 
pinched cylinders with different profiles. This approach 
provides broad possibilities in analyzing the effect of the 
shape of the Fermi surface on observed characteristics of 
layered conductors. 

Let us assume that the function E(p z /po) in the inter- 
val — po < p z < po is described by the expression 



(33) 



where the parameters / and r take values more than unity. 
The model specified by (32) and (33) makes it possible 
to describe a broad class of FSs. 

The Gaussian curvature of the FS is described by the 
expression 





1 


i - 


Ei 


f 




~ 7i 




Po 





K(p z ) = 



m i>z + (P 2 ±/m±)dv z /dp z 



(p]_ + m\v 



2 „2\2 



(34) 



At I = r = 2 the curvature of the FS at its sections by 
the planes p z — and p z = ±p equals: 



k(0) 
k(±Pq) 



6S 1 _ 

u max 

26S 1 

S'min Po 



(35) 
(36) 



where S max and S m i n are the maximum and minimum 
sectional areas of the FS: S maK = 5(0), S m - m = S(±po), 
and SS = S mlix — S m - m = (ir /2)m±r]Vopo. Thus, if the FS 
remains a pinched cylinder (i] ^ 0), its curvature at all 
points of the sections with extremal diameters is finite. 
Similar results are obtained if the tight-binding approxi- 
mation is used to describe the electron energy spectrum. 

For r / 2 and I = 2 the curvature of the FS near 
p z ±po remains finite and k(0) is still described by (35). 
However, the asymptotic behavior of the curvature of the 
FS near p z — ±p is different: 



k(Pz 



-2(r-l) 



5S 1 

Smin Pq 



1-1* 

Pa 



r-2 



(37) 



Thus, for 1 < r < 2 the curvature of the FS has singular- 
ities in these sections. For p z = ±p the curvature n(p z ) 
vanishes at p z = ±po- The corresponding sections of the 




(a) 





- po po 



Pz 



FIG. 2: Profiles of corrugated cylinders described by Eqs.(32) 
and (33), for different values of the parameters r and I; (a) I — 
r = 2; (b) r = 2, 1 > 2; (c) r > 2, 1 > 2; (d) r = 2, 1< I < 2. 



Fermi surface are lines of parabolic points. The larger 
the value of the parameter r, the more the FS near these 
sections resembles a cylinder. 

The anomalies in the curvature of the FS near p z = 
can be described by the model (32) and (33) with r = 2 
and I 7^ 2. Here the curvature of the FS near p z = is 
described by the asymptotic expression 



k{ Pz ) = (/-!) 



SS 1 

Synax Pq 



1-2 



(38) 



For 1 < I < 2 the curvature of the FS has a singularity 
at p z = 0; if I > 2 the Fermi surface near p z = resembles 
a cylinder, and the larger the value of I the closer the 
resemblance. Finally, if r ^ 2 and I ^ 2, we have a 
surface in the form of a pinched cylinder with curvature 
singularities in all the sections with extremal diameters. 
The profiles of the Fermi surface described by (32) and 
(33) are depicted schematically in Fig. 2. 

We shall use our model to obtain the expressions for 
the oscillating function A. The form of the function A 
depends importanty on the characteristic properties of 
the electron energy spectrum. Assuming that FS curva- 
ture on the effective strips is nonzero (r = I = 2) we can 
obtain: 



MO) 




X cos 



X COS ' 



< /2Ci5' max 




heB 


~ 4 


ncS m [ n 




heB 


+ 4 



°° (_-\\n 

(39) 
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Here, So = 27rm_i_£ is the area of the FS cross-section in 
the limit 77 — >■ 0, when the FS turns a pure cylinder; hflo 
is the spin splitting energy. 

In the case when there is an anomaly in the FS curva- 
ture at p z = ±po (( = 2;r ^ 2) the contribution to the 
oscillating function A from the vicinities of correspond- 
ing cross-sections assumes a form: 



AM 



rS \ 1/r r(l/r)™ (-1) 



2j 2 SS 
ncS, 



x cos 



Q 



r e — ' n 1 / 

n=l 



, cos 7TO- „ 
HeB 2r \ Q 



(40) 



where T is the gamma function. Under this condition the 
first term in Eq. (40) retains its form. Correspondingly 
when the local geometry of the FS is characterized with 
the anomalous curvature at p z = (r = 2; I 7^ 2) the 
oscillating function A(0) is described by the expression: 



A(0) 



ISp 
2 7 2 (SS* 



1/1 



T(l/l) ^ (_i)" 



/ ^— ' n 

n— 1 



MO) 



( ncS max tt\ ( tt \ 

x cos —. |— — cos irn—- , (41) 

V h\e\B 21 J V ^ J 

At the same time the expression (40) for A(po) remains 
valid. 

When the effective strip on the FS centered at p z = 
or p z = ±po is close to a cylinder in shape, its contri- 
bution to the function A has the same form as for the 
conductor with two-dimensional energy spectrum. For 
instance, if / 3> 1 the expression (41) can be transformed 
to the form 



00 / 

A (°)^E(- 1 )>«W cos ( 



ncS ms , x \ { tt 
1 — r~, — cos 7m—— 
V h\e\B J \ n 



(42) 

In the case when both parameters (r and I) are large com- 
pared to unity the expression for A contains two terms 
of the form (42). Both terms represent the contributions 
from the quasicylindrical strips on the FS. These oscil- 
lating terms have to differ in period. The difference in 
their periods arises due to the distinction in the extremal 
cross sectional areas. It may be noticeable if the magnetic 
field is not too strong and the inequality ^SS/Sq ^> 1 is 
satisfied. 

In a very strong magnetic field the inequality 
-f 2 5S/So <C 1 is satisfied. Under this condition the dif- 
ference in form between the FS and the cylinder does not 
affect the oscillating function A : 



A 



2^(-l) I V„(6 , )cos(7m 7 2 )cos f tto-^ 

n=l ^ 



(43) 



When a temperature is moderately low (9 ~ 1) we can 
keep only the first term in the sum over "n" in Eqs. (39)- 
(43). Thus, under 9 ~ 1, the dependence of the function 



A on the inverse magnetic field has to be of a harmonic 
type. Because of the factor exp(— 6), the amplitude of 
the oscillations of the function A under 9 ~ 1 is small 
compared to unity even in strong magnetic fields. 

When <§C 1, the form of the oscillations becomes much 
more complicated and their amplitude increases. To eval- 
uate the amplitude of the oscillations under this condi- 
tion one can use the asymptotic formulae following from 
the Euler-Maclaurin summation formula. As a result, 
under this condition we obtain the following estimation 
for the amplitude of the oscillations: 



A max (p ) 
A max (0) 



S 9 \ 1,r 1 



7 2 <5S* 



S Q 9 \ 1/1 1 



7 2 <5S* 



(44) 
(45) 



Here, r and I are the dimensionless parameters whose 
values determine the local geometry of the FS in our 
model (33). In the case when the FS curvature remains 
finite and nonzero at p z — and p z = ±po we have: 



A max (p ) = A max (0) 




MO) 



(46) 



For the same values of the parameters 7 , 5S/ So and 
9 the ratio of the amplitudes of oscillations described 
with the formulae (41) and (39) is of the order of 
("f 2 5S/So9y i ~ 2 ^ 21 . Thus, under typical experimental 
conditions when -f 2 5S/So 3> 1, the oscillations associ- 
ated with a quasicylindrical extremal cross-section of the 
FS have an amplitude much larger than the amplitude 
of the ordinary oscillations described by formula (39). If 
1 < I < 2 one of the principal radii of curvature vanishes 
at all points of the extremal cross-section of the Fermi 
surface at p z = 0. In this case the amplitude of oscil- 
lations described by formula (39) is much smaller than 
the amplitude of the oscillating parts of the contributions 
from the FS ordinary extremal cross-sections. 

When 9 <C 1 the oscillation amplitude may be of the 
order of unity. Charge carriers concentration in layered 
organic metals is far below than in ordinary metals. This 
generates much more favorable conditions for observation 
of the salient features in low temperature quantum oscil- 
lations in thermodynamic quantities. If we take for m± 
and Sq the values obtained in experimental study of the 
Fermi surface of the organic metal (3 — (ET) 2 IBr 2 (Ref. 
H3) (m_L w 4.5m ; S w 1.26 x 10~ 39 g 2 cm 2 / s 2 ; ro is 
free electron mass) the condition (3 /2)1 / 'jy/9 V So/5S ~ 
1 in magnetic fields B - 200 kG and for 6S/S Q w 0.05 will 
be satisfied at temperatures of the order of one Kelvin. 

Now we can return back to our expressions (11) and 
(12) for the elastic moduli and to proceed the analysis of 
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their low temperature anomalies in layered conductors. 
For definiteness, we assume that the FS of the considered 
organic metal is nearly cylindrical in shape near p z = 
and its curvature at p z = ±po is finite and nonzero. 
Under these assumptions the oscillating correction A(0) 
predominates and we can write: 



N 2 _ 7V 2 [l-(l-a (0))X'] 



(1 - a )go 



(47) 



A = Ao 1 + 



1 - 4ttx\\ 



N 2 [l-(l-a (0))K"} 
(1 - a )g 



(48) 

Here, the function K" includes the oscillating term A(0), 
therefore it oscillates itself in a strong magnetic field: 



K" = 



A(0) 



1 + [1 + a (0) - a (0) + A)(0) - 4ttxo7 4 ]A(0) ' 

(49) 

To obtain the expression for K" we used the result (31) 
for the longitudinal part of the magnetic susceptibility 
X|i; the quantity Xo is related to the Landau diamagnetic 
susceptibility. The latter equals — xo/3- 

The sharpest difference between A and Ao is displayed 
at those values of the magnetic field B, at which the 
quantity 1 — 47TX|| goes to zero and the denominator in 
Eq. (49) turns zero at the peaks of the quantum oscil- 
lations. It makes the elastic constant A (or the electron 
contribution to the elastic modulus en) vanish near the 
oscillations maxima. Thus, there can arise lattice insta- 
bility. It is connected with the instability in the magneti- 
zation and arises because of the magnetostriction. Under 
these conditions the Fermi-liquid correlation plays the 
major part. When 47rxo7 4 < 1 the instability can arise 
only under a (0) - a (0) + /3 (0) < 0. 

The lattice instability can, in principle, be displayed in 
the case, when the magnetostriction does not play a part 
(the strain gradient is directed along B), as the function 
K [see Eq. (27)] in the low temperature region can also 
turn to infinity. Such instability was predicted in time 
in Refs. |2^| and However, actually it is probably 
unattainable. The reason is that when A increases the 
quantity l+[l+a (0)-a (0)+/3 (0)-47rK7 4 ] A(0) reaches 
zero before the quantity 1+ [l+a (0) -a (0) +(3 (0)] A(0) . 
Hence, the jump of the magnetic state will take place 
before the constant Ao runs into zero. 



5. CONCLUDING REMARKS 

The point of the present work is that the Fermi-liquid 
correlation among the quasiparticles can strongly affect 
their renormalized DOS. This was analyzed in detail for 
isotropic electron liquid 

0E1E1E3 Here, we gener- 
alize this analysis to cover a broad class of layered con- 
ductors whose FSs are supposed to be axially symmet- 



ric. It is shown as a result of lengthy but straightfor- 
ward calculations that Ni significantly differs from the 
"bare" DOS N{ under the condition that oscillating cor- 
rections A m which arise in the presence of the quantizing 
magnetic field are not small in magnitude compared to 
unity. Under this condition the oscillating part of N£ 
can have singularities in the peaks of quantum oscilla- 
tions. These singularities occur due to the Fermi-liquid 
interaction. The singularities of N£ can cause specific 
features manifested in the oscillations of observables. In 
particular, they can lead to the singularities of the longi- 
tudinal magnetic susceptibility X| at the magnetic fields 
corresponding to the peaks of the quantum oscillations 
of DOS, which leads to the violation of the stability of 
the considered substance (xii < l/47r). A similar effect 
was briefly discussed in Ref. ^lj. Besides we obtain a 
possible singularity in the electron contribution to the 
elastic constant C22 which is connected with the instabil- 
ity in the magnetization and can give rise to the lattice 
instability of a special kind. 

These interesting phenomena can be available for ob- 
servation in organic metals (including substances of the 
a - (BEDT-TTF) 2 MHg(SCH) 4 group) because of the 
specific character of the energy spectra of the charge car- 
riers in these substances. It is too early to draw any con- 
clusions about the local features of the geometry of the 
FS of the majority of layered organic metals, since there 
is a lot to study in the electron energy spectra of such 
materials. It can be assumed, however, that here, as in 
ordinary metals, the FS contains quasicylindrical bands 
or sections with an anomalously large curvature. These 
features of the local geometry of the FS can be created 
(if they are absent) or enhanced by applying an agent 
that changes the shape of the constant-energy surfaces, 
e.g. by applying external pressure along the normal to 
the conducting planes. 

The above analysis shows that the special features in 
the profile of the rippled cylinder, which is the main part 
of the FS of layered organic metals, can substantially in- 
fluence quantum oscillations of thermodynamic observ- 
ables in these materials. The model developed in this 
paper makes it possible to study in detail the observ- 
able manifestations of the local geometry of the FS of 
layered conductors. It resolves some of difficulties that 
emerge when we use the model of tightly bound elec- 
trons. For instance, the characteristic features of the ob- 
servable properties of layered conductors, arising due to 
strong anisotropy of the electrical conductivity can be de- 
scribed and analyzed without passing to the limit rj — » 0, 
which corresponds to a conductor with a two-dimensional 
spectrum of the charge carriers. 

The model specified by (32) and (33) enables us to 
carry out a detailed analysis of quantum oscillations of 
the elastic constants in organic metals and this leads us 
to the conclusion that the local geometry of their Fermi 
surfaces can provide favorable conditions for observation 
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of the above described magnetic and lattice instabilities. 
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APPENDIX 

Here we present in more detail calculation of the renor- 
malized DOS N£. We start from the expansion (23) for 
the matrix elements of the Fermi-liquid kernel. To pro- 
ceed we use the notations: 

vv' 

x Ku'{-Q) S Ta'S Xo;x ' -hcq/eB, (50) 
Iv - fv> 



': = E^<3 |2 "fc»-.(S^FTT 

vv 1 

X Kv'(-q) S <T<T'fix ;x' -hcq/eB, (51) 
vv' 



X $crcr'fixo;x' —hcq/eB 

c' = > 



.90 



E|— |;« |2 "(p»)« |2a,| (p«) 



X ^ff^io^-Scg/eBi (53) 



1 _ 

o' = > 



9o 



EiM^ p|2s| ^ Q|2s ''^ x 



X 5<jo' 5 Xa -,x' a — hcq/eBi 



(54) 



x J;- 2s V2A^+1^ n„/„(q); 

^ - E#^7 p|2s| (^) 

x J;- 2s V27V + 1^ tm„/„(q). 



(55) 



(56) 



The expressions for the renormalized matrix elements 
n* vlJ , (— q) have the form 

where the quantity n* Nlpz {—q) satisfies the relation, fol- 
lowing from the definition (7): 

+ -J2Q 12s1 (pM-2s(^V2NTT)i' s . (57) 
9o s \eB J 

When the Fermi surface is axisymmetric the matrix 
elements n* Nlp ^{— q) do not depend on p z : 

nNi Pz (-q) = Ji VwTT^j . (58) 

The average quantities I s and I' s satisfy a system of linear 
equations which follows from the Eq. (57): 

Is + T,(9ss>Is< + r ss >I' s ,) = N„, 
I's + £(rv sis' +q S s'I' s/ ) = K (59) 

s' 

Using (24) we get the expressions for the coefficients of 
the system (59): 

g ss , = A 2s 5 ss ,+J2^mP l2sl (p m )P l2s ' l ( Pm ) 

m 

x J_ 2s (qR m )J_ 2s ,(qR m ), (60) 

q ss , = B 2s 6 ss ,+Y / ^mQ l2S KPm)Q l2S ' l (Pm) 

m 

X J_ 2s (qRm)J-2s'(qRm), (61) 



TV,' 



= J2 A ' m P l2Sl (Pm)Q l2S \ Pm ) 
m 

X J-2s(QRm)J-2s'(QRm)) 



(62) 



where R m in the radius of the cyclotron orbit correspond- 
ing to the m-th extremal cross section of the FS: 

m± - \ f[ P \*>\( n \V 
) J 

777 I 1 



B 2s — 



27r 2 r 50 



[Q^{p z )} z dp z . (63) 



In obtaining Eqs. (60)-(62) we used the identity concern- 
ing the Bessel functions 



J l-m{x)Ji(x) = 5 l0 . 



(64) 
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When the FS has the unique extremal cross-section (at 
p z = 0) Eq. (64) takes the form: 



1 + A 2s 



+ ^, P\ 2s \Q)J-2 S {qR ex ) Y ^ 



1 + A 



2s 



±+B 2s 



v A , Q l2sl W-2s(qR ex ) x 
s 1 I B 

A Q\ 2s Ko)J-2s%R ex ) Y _ 

l + B 2s 



(65) 



Here: 



s 

Y = J2Q l2Sl (0)J-2s(qRe X )l' s - 



(66) 
(67) 



We can find the quantities X and Y solving the system 
of equations: 



p\2s\ (n\ 

X(l + a q A) + Ya q A> = V N s - U J- 2s (qR ex ), 

1 + A 2s 

^ q\ 2s Uo) 

Xp q A' + Y(l + p g A) = V K?—p. J- 2s (qR ex ). 

1 + B 2s 

(68) 



Here 



a a = 



13, 



A 

+ A 2s 
B* 

J-2s{<lRex) 



^YTt s J - 2s(qRexl 



D 



2 a 



A* 
A 2s 



[pl 2s l(0)] 2 ; B 2 * s -[0 |2s| (0)] 2 . (69) 
As a result we arrive at the expressions: 

u s 1 + A 2s 

Y^K^—U- J- 2s (qR ex ), 



Y 



D 

1 + aA 
D 



Bo 



Bo 



(70) 



+ A 2s 

where the determinant of the system equals: 

D = 1 + (a q + /3 q )A + a q /3 q (A 2 - A' 2 ). (71) 

Sustituting the expressions (70) into the system (65), we 
can calculate the averages: 



h = 



N s A + f3 q (A 2 - A' 2 ) pl 2s l(0) 



l + A 2 



D 



J2 N *' —. X J-2s>(qRe X )- 



12s' 



l + A 2s 
ZV p|2*l(0) 
~~D 1 + A- 



J-2s(qRe X ) 



2s 



X J-2s(qRe X ) V N' s , 7 1 > J_ 2s , (qR ex ), 
, 1 + B 2s < 



(72) 



/: = 



A + a q (A 2 - A' 2 ) Ql 2s l(0) 



l + B 2s 



D 



, 1 



+ A 



l + A 2s 
A' Q\ 2s \(0) 
~D 1 + B 2s 



J-2s(qRe 



pl 2s ' (0) 

J- 2s (qR ex )^2 N s > = — J-2s'(qR ex )- 

, 1 + A 2s i 



(73) 



The asymptotic expressions for the quantities N s and 
N' s also can be found by the formula (24). For the FS 
with the unique extremal cross-section we have: 



N s = - 



• J dp z P a (p z )5 s0 - g AJ (qR ex ) 



X p\ 2s \0)J_ 2s (qR ex ), 
N' s = -g Q A'Jo(qR ex )Q l2sl (0)J-2s(qRe X ). (74) 



Using the Eqs. (57) and (72)-(74), we obtain 

u - u 



50 



1 - Oto 



+ 



[(1 - a ) 2 J 2 (qR ex )][A + [3 q (A 2 - A' 2 )} 



l + (a q + f3 q )A + a q f3 q (A 2 -A' 2 ) 
Here we introduced additional notations: 

A _ Aq 

1 + A ' 



(75) 



a 



l + A 
A =( m \ 



/ P°(Pz)d Pi 



(76) 



Passing to the limit q — > in the Eq. (75) we obtain 



the final result for iV^ . 
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